function [flag,ret_T,edge_T,mid_T]=question_1(alpha,h,range,T_set,v)
%解决问题一

    v=v/60;
    x=0:0.1:range(end);
    q1_time=x/v;
    q1_time=q1_time';

    n=length(q1_time);

    %设置问题一的环境温度
    T0=zeros(size(q1_time));
    for i=1:n
        T0(i)=get_K(T_set,x(i));
    end

    Divide_line=zeros(6,1);%存储各个温区的横坐标位置
    for i=1:4
        m=find(T0==T_set(i));
        Divide_line(i+1)=m(end);
    end
    Divide_line(6)=n;

    u0=25*ones(3,1);
    mid_T=zeros(size(q1_time));
    edge_T=mid_T;
    for i=1:5
        [mid_T(Divide_line(i)+1:Divide_line(i+1)),u0,u(:,Divide_line(i)+1:Divide_line(i+1)),edge_T(Divide_line(i)+1:Divide_line(i+1))]=PDE(alpha(i),h(i),u0,T0(Divide_line(i)+1:Divide_line(i+1)),q1_time(Divide_line(i)+1:Divide_line(i+1)));
    end

    
    [flag]=limit(q1_time,mid_T);%判断是否符合制程要求
    
    figure
    plot(q1_time,mid_T,'red',q1_time,edge_T,'black');
    xlabel('时间/s')
    ylabel('温度/℃')
    legend('中心温度','边界温度')
    
    
    dT=mid_T-edge_T;
    figure
    plot(q1_time,dT,'r')
    xlabel('时间/s')
    ylabel('中心与边缘温度差/℃')
     
    loc=zeros(4,1);
    m=loc;
    loc(1)=111.25;%3中点
    loc(2)=217.75;%6中点
    loc(3)=253.25;%7中点
    loc(4)=304;%8结束
    
    for i=1:4
        [~,m(i)]=min(abs(x-loc(i)));
    end
    ret_T=mid_T(m);
    
end

